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Abstract 



U 
Ph. 

I A model of particle interacting with quantum field is considered. The model 

\^ • includes as particular cases the polaron model and non-relativistic quantum electro- 

CSJ , dynamics. We compute matrix elements of the evolution operator in the stochastic 
approximation and show that depending on the state of the particle one can get the 

■ non-exponential decay with the rate t~2 . In the process of computation a new alge- 

^sj" , bra of commutational relations that can be considered as an operator deformation 

I of quantum Boltzmann commutation relations is used. 

o 

^ ■ 1 Introduction 



For many dissipative systems one has the exponential time decay of correlations. This 
result was established for various models and by using various approximations, see for 
C I ex. IQ. For certain models, in particular for the spin-boson Hamiltonian, also a regime 

with the oscillating behavior was found @] , , 0] . The presence of such a regime is very 
O^' important in the ivestigation of quantum decoherence. The aim of this note is to show 

that for the model of particle interacting with quantum field, in particular for the polaron 
^ I model, one can have not only the standard exponential decay but also the non-exponential 

I decay (as some powers of time) of correlations. 

We investigate the model describing interaction of non-relativistic particle with quan- 
tum field. This model is widely studied in elementary particle physics, solid state physics, 
quantum optics, see for example We consider the simplest case in which mat- 

ter is represented by a single particle, say an electron, whith position and momentum 
q = (gi, q2, qs) andp = (pi,P2,P3) satisfying the commutation relations [qj,Pn] = i^jn- The 
electromagnetic field is described by boson operators a{k) = {ai{k),a2{k),a3{k));a'^{k) = 
(a|(/c), . . . ,al{k)) satisfying the canonical commutation relations [aj{k), al^{k')] = 6jnS{k — 
k'). The Hamiltonian of a free non relativistic atom interacting with a quantum electro- 
magnetic field is 

H = Ho + XHj = J uj{k)a\k)a{k)dk + + Ai^/ (1) 
where A is a small constant, u{k) is the dispersion law of the field, 

Hi = J d^k (g{k)p ■ a\k)e-''"' + g{k)p ■ a{k)e''"') + h.c. (2) 
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Here p ■ a{k) = E?=iPj«i(^), = T^^iP], a^{k)a{k) = E?=i al{k)aj{k), kq = E?=i kjqj. 
For the polaron model the Hamiltonian has the form 

H = J u{k)o){k)a{k)dk + ^p^ + xjd'^k (g{k)a\k)e-'^'i + g{k)a{k)e''"' 

It is different from (||), by a momentum p in the interaction Hamiltonian. For the 
analysis of this paper this difference is not important. 

In the present paper we will use the method for the approximation of the quantum 
mechanical evolution that is called the stochastic limit method, see for example [Q, [^- 
Tl[] . The general idea of the stochastic limit is to make the time rescaling t t/X^ in the 
solution of the Schrodinger equation in interaction picture uj;^^ = e**^oe~**^, associated 
to the Hamiltonian H, i.e. 

with Hj{t) = e^^^°Hje~^^^°. We get the rescaled equation 

and one wants to study the limits, in a topology to be specified, 

ItoC/W (3) 

We will prove that Ut is the solution of the equation 

dtUt = -iHtUt ; f/o = 1 (4) 

The interest of this limit equation is in the fact that many problems become explicitly 
integrable. The stochastic limit of the model (|l|)-(@) has been considered in [|T0|, [|ri|], [|T^, 



HI, [0], E 



After the rescaling t —>■ t/\^ we consider the simultaneous limit A ^ 0, t ^ oo 
under the condition that X^t tends to a constant (interpreted as a new slow scale time). 
This limit captures the main contributions to the dynamics in a regime, of long times 
and small coupling arising from the cumulative effects, on a large time scale, of small 
interactions (A ^ 0). The physical idea is that, looked from the slow time scale of the 
atom, the field looks like a very chaotic object: a quantum white noise, i.e. a (5-correlated 
(in time) quantum field b^j{t, k),bj(t, k) also called a master field. If one introduces the 
dipole approximation the master field is the usual boson Fock white noise. Without the 
dipole approximation the master field is described by a new type of commutation relations 
of the following form |TT 

bj{t, k)pn = {Pn - kn)bj{t, k) (5) 

bj{t, k)bi{t\ k') = 27i6{t - t')6 (u{k) -kp+ 8{k - k')Sjn (6) 

Such quantum white noises can be treated as an operator deformation of quantum Boltz- 
mann commutation relations. Recalling that p is the particle momentum, we see that the 
relation (|^) shows that the particle and the master field are not independent even at a 
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kinematical level. This is what we call entanglement. The relation (Q) is a generalization 
of the algebra of free creation-annihilation operators with commutation relations 

and the corresponding statistics becomes a generalization of the Boltzmannian (or Free) 
statistics. This generalization is due to the fact that the right hand side is not a scalar 
but an operator (a function of the atomic momentum). This means that the relations 

are module commutation relations. For any fixed value p of the atomic momentum 
we get a copy of the free (or Boltzmannian) algebra. Given the relations (|^), (^), the 
statistics of the master field is uniquely determined by the condition 

hj{t,k)^! = 

where \1/ is the vacuum of the master field, via a module generalization of the free Wick 



theorem, see 14 



In Section 2 the dynamically g-deformed commutation relations (|^), (^, ( P^ are 
obtained and the stochastic limit for collective operators is evaluated. In Section 3 the 
stochastic limit of the evolution equation is found. In Section 4 the non-exponential decay 
for vacuum vector in the polaron model is investigated. 



2 Deformed commutation relations 

In this section we reproduce in the brief form the notations and the main results of the 
work ||14 |. 

In order to determine the limit (^ one rewrites the rescaled interaction Hamiltonian 
in terms of some rescaled fields axj{t, k): 



J Hi (^^) = J d^kp{g{k)ax{t, k) + g{k)a{{t, k)) + h.c. 



where 



ay,^j{t,k) := ^e'^^Oe'%j(fc)e-'^^° = le"*^ ("('=)-^^+^^'')e''=%-(A;) 
A A 

It is now easy to prove that operators a\j(t, k) satisfy the following g-deformed module 
relations, 

ax,j{t,k)a{ {t',k') = 



i{jt', k')ax,j{t, k) -qxit- t', kk') + ^qx(t- t', uj{k) - kp + 6{k - k%n (7) 

axj{t, k)pn = {Pn - kn)axj{t, k) (8) 

where 

g,(t-t',x) = e-^'#^ (9) 



is an oscillating exponent. This shows that the module g-deformation of the commutation 
relations arise here as a result of the dynamics and are not put artificially ah initio. For a 
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discussion of g-deformed commutation relations see for example |16|. Now let us suppose 
that the master field 

bj{t,k) = \imax,jit,k) (10) 

A — >0 

exist. Then it is natural to conjecture that its algebra shall be obtained as the stochastic 
limit (A — > 0) of the algebra (0), (|]). Notice that the factor qx(t — t',x) is an oscillating 
exponent and one easily sees that 

lim qx{t, x) = , lim ^ qx{t, x) = 27i6{t)6{x) (11) 

A — >0 A — >0 /\ 

Thus it is natural to expect that the limit of (^ is 

bj{t, k)pn = {Pn - kn)bj{t, k) (12) 
and the limit of (|^ gives the module free relation 

bjit, k)bl{t', k') = 2Tx5{t - t')5 [uj{k) -kp+ 6{k - k')5jn (13) 

Operators axj(t, k) also obey the relation 

ax,j{t, k)ax,n{t', k') = ax,n(t', k')ax,j(t, k)ql^{t - t' , kk') (14) 
In what follows we will not write indexes j, n explicitly. The relation (|T^) should disappear 



after the limit, see [|Tj]. In fact, if the relation (|Tj) would survive in the limit then, 
because of (p!l|), it should give 6(t, k)b(t', k') = 0, hence also b^{t, k)b^{t', k') = 0, so all the 
n-particle vectors with n > 2 would be zero. 



3 Evolution equation 

Let us find stochastic differential equation for the model we consider. In the introduction 
we claimed that the stochastic limit for the Shrodinger equation in interaction picture will 
have the form (^: dtUt = —iHtUt- But in this equation both Ht and Ut are distributions. 
We need to regularize this product of distributions. In the present section we will make 
the following regularization: roughly speaking we replace Ht by Ht+o + const. 

We investigate the evolution operator in interaction picture Ut^\ We start with the 
equation 



+ 



ft+dt j-ti \ 



(A) 

t 



where rft > 0. We get for d^/f ^ = u[%^ - f/f ^ 

rff/f ^ = {[-iX] J'^'' Hj{t,)dh + i-zXy f^'"'' dt, f^' dt2Hj{t,)Hj{t2) + . . .) Uj; 
Let us make the rescaling t t/X"^ in this perturbation theory series. We get 

=((-.) if, (^) + 



(A) 



4 



To find the stochastic differential equation we need to collect all the terms of order dt 
in the perturbation theory series (0). Terms of order dt are contained only in the first 
two terms of these series. Let us investigate the first two terms. For the first term of the 
perturbation theory we get 

dh-Hi (^^j = dh j dk <(g{k){2p + k)ax{ti, k) + g{k)a{{U, k){2p + k)) 

(16) 

In the stochastic limit A ^ this term gives us 

dk (g{k){2p + k)dB{t, k) + g{k)dB\t, k){2p + k 

where the stochastic differential dB{t, k) is the stochastic limit of the field ax{t., k) in the 
time interval (t, t + dt): 

rt+dt rt+dt 

dB{t,k) = \im / dTax(T,k)= / dTb{T,k) 

We will prove that the stochastic differential dB{t, k) and the evolution operator Ut are free 
independent. In the bosonic case independence would result in the relation [dB(t, k), Ut] = 
0. From this relation follows that {X dB{t, k)Ut) = for arbitrary observable X. In the 
case of Boltzmannian statistics we get the same relation: the (free) independence means 
that roughly speaking dB{t, k) kills all creations in Uf. We have the following 

Lemma. The stochastic differental dB{t, k) and the evolution operator Ut are free 
independent. This means that for an arbitrary observable X 

{XdB{t,k)Ut) = VX 

Here (■) is the stochastic limit of the vacuum expectation of boson field (that acts as a 
conditional expectation on momentum of quantum particle p). 

We will prove this result by analizing of the perturbation theory series. We have 



rt+dt / /■* 1 / f 1 

{X dB{t, k)Ut) = lim(XA draxir, k) [l + (-z) dt^ - Hj i^^ 



+ ...+ 



+Hr/;*,.../;"-di„ii/,(^)... iff, (I ) + ...)) 



where jHj (^j^^ is given by the formula (|16|). Here limA^o-^A = X. Let us analize the 
A^-th term of perturbation theory. The iV-th term of perturbation theory is the linear 
combination of the following terms (we omit integration over k, kn) 

rt+dt rt ptN-l 

{Xx dTax{T,k) dti... dt^ a^'xiti.ki) . . .a"^ {tN.kN)) 

Jt Jo Jo 

Let us shift ax{T,k) to the right using dynamically g-deformed relations. In the follow- 
ing we will use notions of the work |]T^. Let us enumerate annihilators in the prod- 
uct a^^ {ti, ki) . . . a'l^ {tN, kj^) as ax{tm,km), j = I,... J, and enumerate creators as 
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a\{tm'., km'.), j = 1, . . . I, I + J = N. This means that if = then a'j^ltm, km) = 
axitm , km ) for m = rrij (and the analogous condition for Em = !)• 

We will use the following recurrent relation for correlator (analogous formula was 
proved in the work [p!4| ): 

/■t+dt rt /"tjv-l 

{Xx dr / dti... dtN ax{T,k)al'^{ti,ki) . . .al'^{tN,kN)) = 



^0 



^ rt+dt rt rtjv-1 , 

\S{k-km'.){Xx dr dti... dtjy al^{ti, ki) . . .al{tm'., km'.) ■ ■ ■ al'^ {t^, k^)) 

^ Jt Jo Jo J J 



[t - tm'. , Uj{k) -kp+ i/c^^ Yi ?A ^ - tm'j , kkm,) Yl - , kkm'. 



mi>m'. 



n Qx^i^ -^m,,kkm,) n Qxir - tm'^, kkm'J (17) 



mi<m'. 



Here the notion d\ means that we omit the operator a\ in this product. 
The right hand side of the equation (|T^) is equal to 

y^5{k-km'){Xx / dti... dtNaf{ti,ki)...d\{tm'.,km'.)-.-al'^{tN,kN)) 
~i ^ Jo Jo J J 

^qx (^-tm'^,uj{k) -kp+ ^k^^ Yl q^^ [-tm'^,kkm,) Yl ^A (^-tm'pkkm'^ 



mi>m' m'.>m'. 



n (lx^i.-im,,kkm,) n Qxi-tm'^, kkm'J 
mi<m'. m'^<m'. 

rt+dt / 1 \ / 

J dr qx [t, uj{k) - kp + -k'^j Yl Qx^ kkm,) Yl kkn 
n Qx^i'^^kkmJ Yl (lx{r,kkm') 



m'. <m'. 
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(we use that qx is an exponent). The first three lines of this formula do not depend on 
r and the last two lines do not depend on ti, . . . ^t^. Therefore the stochastic limits for 
these values can be made independently (the limit of product is equal to the product of 
limits). It is easy to see that the stochastic limit for the multiplier that depends on r (of 
the last two lines) is equal to zero. This finishes the proof of the Lemma. 

The second term of the perturbation theory series is equal (up to terms of order {dt)"^) 

= j'J" dt, j[" dh j d.k\g{kn2p + k fl-^ e-^ ("(^^-P+i -.") 

due to g-module relations on ax(t, k), p. Performing integration over ti, t2 and using the 
formulas 

rt+dt rti 1 .ti-io rt+dt rO 

dti / dt2TT7e-*^"= / dti / c^re^™ 
Jt Jt J-ti/x^ 



f° dte'''' = — = 7i5{x) - i P.P.- 

J — oo X — tO X 



we get for the second term 

-idt J dk\g{k)\\2p + 

Let us denote 



uj{k) -kp+^k'^-iO 



{9\9)-{p) = -^ I dk\g{k)\\2p + ky ^ 



uj{k) -kp + ^k^-iO 



Combining all the terms of order dt we get the following result: 

Theorem. The stochastic differential equation for Ut = lim f///'^_^2 have a form 

dUt = (^-i J dk(g{k){2p + k)dB{t,k) + g{k)dB\t,k){2p + k)) - dt {g\g)4p)jUt (18) 



The equation (|T8D can be rewritten in the language of distributions as 
^ = (^-tJdk(gik)i^P + mt,k)+gik)b\t,k)i2p + k))- (19) 

Here we uderstand the singular product of distributions b(t, k)Ut in the sense that (|19D is 
equivalent to (|18]). We have to stress that dB{t,k) = J^'^'^^ dTb{r,k) ^ b{t,k)dt and we 
can not obtain ([l^) dividing ([18|) by dt. 

4 Non-exponential decay 

Let us investigate the behavior of {Ut) using stochastic differential equation (|^). We get 

{dUt) = (f(-z) / dkg{k)i2p + k)dB{t,k)-dt{g\g)4p)]Ut) 



Using the free independence of dB(t, k) and Ut we get 

j^{Ut) = {j^Ut) = -{g\g)-{p){Ut) 
Because ?7o = 1, we have the solution 

{Ut) = e-*^^!^)-^^') 

In this section we calculate the matrix element {X\Ut\X) where X = /(p) $ in the 
momentum representation and $ is the vacuum vector for the master field. This matrix 
element is equal to 

{X\Ut\X) = |dp|/(p)|V*(^l5)-(^) (20) 



We investigate the polaron model when u!{k) = 1. For this choice of uj{k) we get 

u;{k) -kp + ^e = l- + _ p)2 
One can expect non-exponential relaxation when supp 

f{p) C {bl < V2}. In this case 
Iie{g\g)^{p) = and there is no dumping. All decay in this case is due to interferention. 

We will use the approximation diam supp g{k) ^ diam supp f{p)- Physically this 
means that the particle is more localized in momentum representation than the field. This 
assumption seems natural because the field's degrees of freedom are fast and the particles 
one are slow. Under this assumption we can estimate the matrix element (pOl). We will 
prove that in this case there will be polynomial decay. 

For \p\ < V2 we get 

{g\g)-{p) = -^ fdk\g{k)\\2p + kf 



I - Ip"" + lik - py 

= -2i J dk\g{k)\'-t{h + h); 



6 J dk\g{k)\'^p{k — p)- 



2^ 

1 



Here only Ji and I2 depend on p and therefore can interfere. Let us find the asymptotics 
of {g\g)~{p) on p (we investigate the case when p is a small parameter). 

We will use the following assumption on g{k): let g{k) be a very smooth function. 
This means that \g{k)\'^ = XF{Xk), F{k) > is compactly supported smooth function, A 
is a small parameter. Let us consider the Taylor expansion on the small parameter p 

XF{Xk) = XFiXik - p)) + X^Y^Pi^PiKk -P)) + .... 

We get that XF{X{k — p)) is a leading term with respect to A. Taking A ^ we get 
that we can use \g{k — p)P instead of \g{k)\'^ in the formulas for Ji and I2 for sufficiently 
smooth g{k). Let us calculate Ji a d I2. Using assumptions considered above we get 



h = (-2 + Wp') f dk \g{k-p)\' 2 ^ 



+ h{k — pY 



(-2 + W) J dk \9m'Y^2 -P'jdk \g{k)\'—^^ 

2 I X ~|~ ~ 



(1 + 



l2=pQ, Q = Q j dk\g{k)\^k- 



2'^ 



We get 



i9\9)-{p) = -2^ / dk \gik)\' + 2ijdk | |2-_L_ _ _ 
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A = 10 fdk \g{k)\'-^ - fdk \g{k)\'- ^ (21) 



We get for X(t) = {X\Ut\X) 

x{t)= J dp|/(p)|V*(^i^)-(^) = 



Let us estimate this integral for f{p) = e , B » 1. Let us consider for simplicity the 
case Q = (for example g{k) is spherically symmetric). In this case the integral is equal 
to 



3 

io \B-iAtJ 



We get that for large t the decay of the matrix element X{t) = {X\Ut\X) is proportional 

3 

to {At)~^ where A is the functional of the cut-off function given by (pl|). 

To conclude, in this paper we obtain that in the polaron model for some (symmetric 
and very smooth) cut-off functions we have the polynomial relaxation, the matrix element 
being proportional to t~2. The dependence on the parameter B that corresponds to the 
size of the support of the smearing function f{p) of quantum particle in the momentum 
space for large t is not important. We can say that particles with large momentum decay 

3 

exponentially and the particles with small momentum decay as t~2 and the decay for 
large t does not depend on the smearing function /(p). 
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